Search for: All records

Understanding the connections between galaxy stellar mass, star formation rate, and dark matter halo mass represents a key goal of the theory of galaxy formation. Cosmological simulations that include hydrodynamics, physical treatments of star formation, feedback from supernovae, and the radiative transfer of ionizing photons can capture the processes relevant for establishing these connections. The complexity of these physics can prove difficult to disentangle and obfuscate how mass-dependent trends in the galaxy population originate. Here, we train a machine-learning method called Explainable Boosting Machines (EBMs) to infer how the stellar mass and star formation rate of nearly 6 million galaxies simulated by the Cosmic Reionization on Computers project depend on the physical properties of halo mass, the peak circular velocity of the galaxy during its formation historyv_peak, cosmic environment, and redshift. The resulting EBM models reveal the relative importance of these properties in setting galaxy stellar mass and star formation rate, withv_peakproviding the most dominant contribution. Environmental properties provide substantial improvements for modeling the stellar mass and star formation rate in only ≲10% of the simulated galaxies. We also provide alternative formulations of EBM models that enable low-resolution simulations, which cannot track the interior structure of dark matter halos, to predict the stellar mass and star formation rate of galaxies computed by high-resolution simulations with detailed baryonic physics.

The analytic galactic wind model derived by Chevalier and Clegg in 1985 (CC85) assumes uniform energy and mass-injection within the starburst galaxy nucleus. However, the structure of nuclear star clusters, bulges, and star-forming knots are non-uniform. We generalize to cases with spherically-symmetric energy/mass injection that scale as r−Δ within the starburst volume R, providing solutions for Δ = 0, 1/2, 1, 3/2, and 2. In marked contrast with the CC85 model (Δ = 0), which predicts zero velocity at the centre, for a singular isothermal sphere profile (Δ = 2), we find that the flow maintains a constant Mach number of $\mathcal {M}=\sqrt{3/5} \simeq 0.77$ throughout the volume. The fast interior flow can be written as $v_{r \lt R} = (\dot{E}_T/3\dot{M}_T)^{1/2} \simeq 0.41 \, v_\infty$, where v∞ is the asymptotic velocity, and $\dot{E}_T$ and $\dot{M}_T$ are the total energy and mass injection rates. For $v_\infty \simeq 2000 \, \mathrm{km \, s^{-1}}$, $v_{r\lt R} \simeq 820 \, \mathrm{km\, s^{-1}}$ throughout the wind-driving region. The temperature and density profiles of the non-uniform models may be important for interpreting spatially-resolved maps of starburst nuclei. We compute velocity resolved spectra to contrast the Δ = 0 (CC85) and Δ = 2 models. Next generation X-ray space telescopes such as XRISM may assess these kinematic predictions.